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We show that density- weighted moments of the dissipation rate, ej, averaged over a scale in supersonic 
turbulence can be successfully explained by the She and Leveque model [Phys. Rev. Lett. 72, 336 (1994)]. 
A general method is developed to measure the two parameters of the model, 7 and d, based directly on their 
physical interpretations as the scaling exponent of the dissipation rate in the most intermittent structures (7) 
and the dimension of the structures (d). We find that the best-fit parameters (7 = 0.71 and d = 1.90) derived 
from the e; scalings in a simulation of supersonic turbulence at Mach 6 agree with their direct measurements, 
confirming the validity of the model in supersonic turbulence. 

PACS numbers: 



Supersonic turbulence is ubiquitous in the cold interstellar 
medium fll] and is believed to play a crucial role in the pro- 
cess of star formation 10, [3]. If supersonic turbulence were 
characterized by universal statistics, as often assumed for in- 
compressible turbulence, the universality could constitute the 
foundations for a statistical theory of star formation. In this 
Letter we focus on the statistical properties of the most inter- 
mittent structures (MISs) of supersonic turbulence. 

The theory of fully developed turbulence assumes that the 
scaling behavior of small-scale fluctuations in the inertial 
range is flow-independent, e.g., the moments of the velocity 
difference, {Sv{1)p) = {{v{x + I) — v{x))p) oc Z^^, have uni- 
versal scaling exponents, (p. The universal state for fully de- 
veloped incompressible turbulence proposed by Kolmogorov 
in 1941 (K41) [4], with (p = p/3, has been shown to deviate 
significantly from (p measured in both experiments and nu- 
merical simulations, at p > 3. This discrepancy is due to spa- 
tial fluctuations in the dissipation rate, neglected in the K41 
theory |6]. The scaling exponents of the average energy dis- 
sipation, over a scale I (see eq. (3)) give corrections to the 
K41 theory, referred to as intermittency corrections flQ]- A 
careful study of fluctuations in the energy dissipation is essen- 
tial for understanding intermittency in turbulence. 

The intermittency model by She and Leveque ^ (hereafter 
the SL model) is based on a hierarchy of dissipative struc- 
tures of different intensity levels, characterized by the ratios, 

e[^'' — {ef^^)/{ef), of successive moments of e;. With in- 

(p) 

creasing order, p, ef represents structures of increasing in- 
tensity and ef°^ corresponds to the MISs. By invoking a 
hypothetical "hidden symmetry" that relates this hierarchy of 
structures to the most intermittent ones, the model predicts 
the scaling exponents, r^, of the energy dissipation moments, 
(ef ) oc V'', of all orders, p. 



-7P + 7(1-/3P)/(1-/?). 



(1) 



The parameter 7 is the scaling exponent of the dissipation rate 
in the MISs, ef"'' oc limp^oo Up+^^'^p oc l^'' , and (3 is re- 
lated to 7 and to the Hausdorff dimension, d, of the MISs by 
7/(1 — /?) = D — d, where D ^ i for three-dimensional (3D) 
turbulence. The physical meaning of this relation will be ex- 



plained later. This model is very successful in predicting ^p in 
incompressible turbulence with high accuracy. 

In this Letter, we study the fluctuations of the dissipation 
rate in supersonic hydrodynamic (HD) turbulence using nu- 
merical simulations. We show that the simulation results for 
the scaling exponents, Tp, are well represented by eq. (1), 
suggesting the SL formulation for the scaling behavior of the 
dissipation rate, originally proposed for incompressible tur- 
bulence, may be applied to supersonic turbulence as well. We 
present a method to directly measure the parameters 7 and 
d according to their physical interpretation, which is general 
and not limited to the supersonic regime of interest here. This 
method can be used to test the validity of the physical inter- 
pretation of the SL model in any of its applications. For super- 
sonic turbulence, we find that the parameters derived with this 
method are in excellent agreement with the values that best fit 
Tp, which confirms the physical interpretation of the model. 

Instead of directly investigating the statistics of the dissipa- 
tion rate, most studies of this model are primarily concerned 
with the scaling exponent, ^p, of the velocity difference, 6v{l). 
The model predicts Cp (1 - 7)p/3 + 7(1 - I3^'^)I{1 - /?), 
which follows from the refined similarity hypothesis, 5v{l) ~ 

and eq (1). Assuming that the largest available ki- 
netic energy in the strongest structures is ~ J7^, with U be- 
ing the rms velocity, and that the timescale in these struc- 
tures follows the usual Kolmogorov scaling, <; oc P/^, She 
and Leveque argued that e^°°^ C/^/i; oc l^"^/^, i.e., 7 = 
2/3 i2lll . With this 7 and with d = 1, corresponding to fil- 
amentary dissipative structures, the values of C^p predicted by 
this model agree with experimental results of incompressible 
turbulence with an accuracy of about 1 % Ist] . 

Although not directly measured from the MISs, (cf. 
7 = 2/3 has been adopted in most applications of the model to 
incompressible ifioll and supersonic MHD turbulence flll[l2ll . 
The dimension d was obtained either from the assumed geom- 
etry of the MISs (d = 2 for current sheets or shocks in MHD 
or supersonic turbulence), or from the best fit to the numerical 
velocity structure functions Il2lll3[l . These works have shown 
that the SL model with 2D MISs is generally consistent with 
simulations of MHD and highly compressible turbulence. 

However, a strict verification of the validity of this model 
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requires a demonstration that the parameters that fit Tp also 
have the declared physical meaning, otherwise the agreement 
between the model and the simulations may be a mere coin- 
cidence. In the present work we thus obtain 7 and d both by 
a direct measurement from their interpretations, and by fitting 
Tp. We are also interested in deriving the dimension of the 
MISs in supersonic HD turbulence where about 1/3 kinetic 
energy dissipates in dilatational modes and 2/3 in solenoidal 
modes, with strongest shocks generally coinciding with the 
locations of strongest shear and vortices. 

Measuring r^.-We take the 1024'^ simulation of supersonic 
HD turbulence for isothermal ideal gas with a rms Mach num- 
ber of 6 from reference Il4ll . The simulation employs the 
piecewise parabolic method to solve the Euler equation fi3\ . 
We focus on the statistics of e/. The dissipation rate per unit 
mass at a given position and time is calculated by jisll . 

e(x, t) = i2Re)-\d,v, + d,v, - {2/i)5,jdkVkf (2) 

where Re is the effective Reynolds number controlled by nu- 
merical dissipation. We compute the velocity gradients at the 
resolution scale and assume Re is constant. 

We calculate the average dissipation rate, e/(x, t), at a scale 
I around x, from the definition given in Jsl] (generalized to 
account for density fluctuations), 

e;(x,t) = ^^^^ \)V{1) j ^(^ + ^''*)^(^ + ^''^)'^^' (3) 

|x'|<! 

where V{1) — Anl^ /"i is the volume of a spherical region of 
size Z, and pi{x,t) = /3(x + x',t)dx' is the 

average density of that region. For convenience, we divide 
the simulation box into cubes (instead of spheres) of different 
sizes in our computations. 

The moments, (e[), of e;, can be evaluated by 

= ^ / ^f(x,i)pKx,i)dx (4) 

where V is the total volume of the system and p is the overall 
average density. We have used the Favre fl^ density weight- 
ing factor, pi/p, to account for the density variations in com- 
pressible turbulence. With this density weighting, the first or- 
der moment, i.e., (e;), is independent of scale I, as follows 
from eqs. (3) and (4), resulting in ti — (see Fig. 1). This 
suggests that, if the refined similarity hypothesis applies to su- 
personic turbulence, the density-weighted third-order veloc- 
ity structure function in compressible turbulence would have 
<^3 = 1, an exact result for incompressible turbulence, known 
as Kolmogorov's 4/5 law |14]. Since the SL formula, eq. 
(1), gives Ti = 0, it is appropriate to compute the density- 
weighted moments from the simulation data and compare with 
the model. The Hausdorff dimension of the MISs we obtain 
from fitting Tp is thus density-weighted in the sense of density- 
weighting in eq. (4). For a valid comparison, we will include 
density-weighting in our direct measurement of d. 
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FIG. 1: Scaling exponents, Tp, of the dissipation rate, ej, averaged 
over 9 snapshots. Error bars indicate snapshot-to-snapshot varia- 
tions. The best fit gives 7 = 0.71 and d = 1.90. 

We calculate moments of e; from eq. (4) for 9 snapshots of 
our simulation, covering more than 5 dynamical times. We ob- 
tain Tp from least-square fits to the log]^Q((ef })-log]^o(0 curves 
in each snapshot. The results are shown in Fig. 1, where the 
data points and the error bars are, respectively, the average 
exponents and the standard deviations over the 9 snapshots. 
The eiTor bars are negligible for p < 2, meaning that there 
are little snapshot-to-snapshot variations for the exponents at 
low orders. The scatter increases with the order, and the error 
bar at p 4 is already significant (7%). We find that, start- 
ing fromp = 4, the log]^Q((e^})-log]^Q(Z) curves are no longer 
well fit by straight lines and thus we only show results up to 
the 4th order Note, however, that the 4th order moment of e; 
coiTesponds to 12th order moment of 6v{l). 

Comparing with eq. (1), we find that the SL model with 
7 = 0.71 and /3 — 0.35 (d = 1.90) gives an excellent fit to 
the numerical data. The fit shows that the SL model can be 
successfully applied to the density-weighted statistics of the 
dissipation rate in supersonic turbulence. As mentioned ear- 
lier, a demonstration that the best-fit parameters indeed carry 
their physical meaning is needed to verify the validity of the 
model. To this end, we directly measure the parameters from 
the simulation data. A fairly large range of parameter pairs, 
(0.67-078) for 7, and a corresponding range of (2.04-1.60) for 
d, can give acceptable (but poorer) fits to the numerical re- 
sults for Tp within the 2 cr eiTor bars. If the SL model works, 
a direct measurement would fix these parameters. 

Measuring 7.-We obtain 7 directly by measuring the aver- 
age dissipation rate profile around the MISs. We first locate 
the MISs in the computational domain, by selecting cells with 
dissipation rate larger than a given threshold, eth, set to be 
close to the maximum dissipation rate over the domain, e„i. 
We will call these cells the dissipation peaks. We then use cu- 
bic boxes of different sizes, I, to cover each peak, and evaluate 
the average dissipation rate in each box, ep{l), through eq. (3) 
(with x at the peak). Taking the average over all peaks, we 
obtain an average dissipation rate, {ep{l)), as a function of the 
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FIG. 2: Average dissipation rate, {ep{l)), (normalized to average 
maximum, {em)) around dissipation peaks selected by thresholds, 
1/3, 1/2 and 2/3 em- Best-fit lines give 7 = 0.705, 0.710 and 0.710, 
respectively. The inset shows convergence of measured 7 with eth - 



cube size, I. The slope of the profile is expected to be the ex- 
ponent, 7. We calculated the average ep{l) in two different 
ways: with and without the average density Pp{l) in a box of 
size I as a weighting factor. We find little difference between 
the slopes obtained in the two ways, implying a weak corre- 
lation between the dissipation rate and the density around the 
peaks We increase the threshold and check whether the 
slope of ep{l) converges. The converged slope is the parame- 
ter 7 that we pursue. 

Our result is shown in Fig. 2 for three thresholds, 1/3, 1/2 
and 2/3 e„i- The three curves are the average over all the peaks 
in the same 9 snapshots used to calculate Tp. The profiles are 
approximated well by power laws (except atl — 1, i.e., at the 
resolution scale) and we find 7 — 0.705, 0.710 and 0.710, re- 
spectively, for the three thresholds. This value is very close 
to 2/3 proposed by She and Leveque, suggesting that the Kol- 
mogorov scaling for the timescale in the MISs, ti oc 
applies also to supersonic turbulence. The measured 7 con- 
verges to 0.71 at the threshold of 1/2 e^, which concides with 
the value obtained from the best fit to Tp. Besides showing 
the applicability of the SL model to supersonic turbulence, 
we have thus verified that the parameter 7 carries the precise 
physical meaning in the model. 

Measuring d-The Hausdorff dimension, d, enters the SL 
model through the following argument, which also provides 
an explanation of the relation 7/(1-/3) = D — d. In the 
limit p — *■ cxD, the contribution to (e|'} would be primarily 
from the MISs at scale I. Since the average dissipation rate 
in regions of size I containing MISs goes like I ^ , we have 
(ef ) oc l^'''PP{l), in the limit p ^ 00, where P{1) is the 
probability of finding a region of size I that hosts a dissipative 
structure of the highest level at scale I. A geometric consid- 
eration suggests that P{1) cx l^^"^ if the dimension of the 
MISs is d [18]. This gives (ef) cx l-iv+D-d p ^ 00. It 
immediately follows from eq. (1) that 7/(1 — (3) = D — d. 

Directly measuring the Hausdorff dimension of the MISs, 
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FIG. 3: Probability of finding cubes of size I with average dissipation 
rate larger than 1/6, 1/4 and 1/3 emZ~"'^^. Scaling exponents of 
P{1) correspond to d = 1.92, 1.90 and 1.90, respectively. The inset 
shows convergence of measured d with eth- 



e.g., using a box-counting method, is challenging. Here we 
take a simpler approach: we compute the probability, P{1), 
of finding an MIS in a cube of size I, and derive d from the 
scaling of P{1) with I, based on the physical argument given 
above We need a criterion to judge whether a cube of size I 
in the simulation box contains an MIS at that scale. Based 
on the log-Poisson version of the SL model jsi [l9|] . we find 
that the appropriate criterion is that the cube in question has 
an average dissipation rate larger than a threshold that scales 
like £thl~'' with L The factor l"'^ accounts for the decrease of 
the average dissipation rate in the MISs with scale. We will 
let eth approach in the simulation box. 

The chosen threshold is justified as follows. The SL 
model is equivalent to a log-Poisson distribution for e/, i.e., 
P{ei)dei - {l/L)^D-^^j:^^,^PL{h-i{eiP) - HeL^) - 
n\n{P))d\n(ei) where L is the integral scale, e the overall 
dissipation rate, and X — (D — d)\n{L/l) jsl I20I1 . The dis- 
tribution, Pi(ln(ei/e)), of the dissipation rate, e^, at L, de- 
pends on the driving force and is thus non-universal. At the 
integral scale, is approximately equal to e, thus Pl{x) is 
supposed to be narrow and decrease very rapidly with in- 
creasing X- If the SL model is valid for supersonic turbu- 
lence, the probability, P{1), of finding a region of size I with 
ei > ethl ' is given by the cumulative probability P{1) = 
(Vi)(^-'')S;f=o^/in(.,./e-)-.in(L)-„i„(,)^i(^)'^^ where 
I and L are in units of the resolution scale. Due to the rapid de- 
crease of Pl{x) with X, for a large eth the contribution from 
the n-th term to P(l) decreases quickly with n because the 
lower integral limit increases with n (since (3 < 1). As eth in- 
creases, the contribution would be more and more dominated 
by the n = term, which is oc 

liD-d)^ Therefore, in our 
measurement, we increase eth and check whether the scaling 
exponent of P{1) converges. If the SL model is coiTect, the 
converged exponent is expected to be the codimension D — d 
of the MISs and to agree with that derived from Tp. To be con- 
sistent with the density-weighting in the dissipation rate mo- 
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ments, each cube that satisfies the criterion is given a weight- 
ing factor proportional to the average density in the cube. 

Fig. 3 shows the scaling of P{1) with I (averaged from the 9 
snapshots), for 3 different e^^: 1 /6, 1/4, and 1 /3 £„. For the 3 
choices of cth, the scaling exponents of P{1) are, respectively, 
1.08, 1.10, and 1.10 (meaning d = 1.92, 1.90 and 1.90). As 
eth increases, d converges to 1.90, which again agrees exactly 
with the best-fit value from the scahng of the dissipation rate 
moments. Together with the measurement of 7, this estimate 
of d validates the extension of the SL model to supersonic 
turbulence, and confirms the validity of the physical interpre- 
tation of the parameters. 

In conclusion, we have studied the statistics of energy dis- 
sipation in supersonic HD turbulence at Mach number M = 6 
using numerical simulations. We have computed the scaling 
exponents, Tp, of density-weighted moments of the dissipa- 
tion rate, ei, averaged over a scale and found that the SL 
intermittency model with 7 = 0.71 and d = 1.9 gives an 
excellent fit to the measured Tp. We stress that, with den- 
sity weighting, n = 0, suggesting a linear scaling for the 
density-weighted 3rd order velocity structure function. We 
have developed a general method to directly measure 7 and d, 
which provides a validity test of the physical interpretation of 
the model. We have shown that the parameters measured di- 
rectly are exactly equal to the values that best fit Tp. We have 
thus verified that the SL model can be successfully applied to 
supersonic turbulence. Investigations with other Mach num- 
bers, especially larger ones, would advance our understanding 
of the energy dissipation in supersonic turbulence. At large 
enough M, there may exist an asymptotic state (possibly al- 
ready reached at M = 6), where the scaling of the energy 
dissipation rate (i.e., Tp) would be universal and independent 
of M, and so would 7 and d. This conjecture is based on the 
observation that, at M > 6, an equilibrium in kinetic energy 
partition between the solenoidal modes (2/3) and the potential 
modes (1/3) is always established (for an isothermal equation 
of state), regardless of their energy ratio in the driving force. 
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